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Abstract 



The analysis of meson correlations by Hanbury-Brown-Twiss interferom- 
etry is tested with a simple model of meson production by resonance decay. 
We derive conditions which should be satisfied in order to relate the measured 
momentum correlation to the classical source size. The Bose correlation ef- 
fects are apparent in both the ratio of meson pairs to singles and in the ratio 
of like to unlike pairs. With our parameter values, we find that the single 
particle distribution is too distorted by the correlation to allow a straight- 
forward analysis using pair correlation normalized by the singles rates. An 
analysis comparing symmetrized to unsymmetrized pairs is more robust, but 
nonclassical off-shell effects are important at realistic temperatures. 
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I. INTRODUCTION 



An important observable for studying high energy collisions is the momentum correlation 
between emitted mesons. In the statistical limit the correlation function depends on the 
Bose symmetry of the particles as in the Hanbury-Brown-Twiss (HBT) analysis of photon 
correlations from stars (for a review see [J]). Under the statistical assumptions, the meson 
correlation directly measures the size of the meson interaction zone, and thus indirectly 
reflects the hadronic dynamics prior to the meson freeze-out. In particular, the correlation 
between pions and between kaons has been analyzed to give quantitative data on the mean 
square source dimensions [0[J. 

However, the analysis is dependent on the statistical assumptions, raising the question 
of their validity under realistic conditions of meson production. In this work we develop 
a simple model of meson production to test the HBT analysis. Within this model we can 
also examine the concept of a classical source. Before defining the model, we shall briefly 
review the HBT analysis to establish the notation. The basic observables are the one- and 
two-particle differential cross sections, da^/d 3 q and da^ / 'd 3 qid 3 q 2 , within some class of 
events. Writing the cross section for these events as o"o, we define the doubles-to-singles 
correlation function as |7j 



The three- vectors q l and q 2 denote the three-momenta of the two mesons. To facilitate later 
on the interpretation of the correlation function, we chose the arguments to be the average 
and relative momentum, q av = (q 1 + q 2 )/2 and q Icl = q 1 — q 2 , respectively. 

In the theoretical model, it will be more convenient to express the correlation function 
using decay rates rather than cross sections. The rates are related to cross sections by a flux 
factor. In terms of the total decay rate W and the differential decay rates W ( - 2 \q 1 , q 2 ) = 
dW/d 3 q\d 3 q 2 and W^'{q) = dW/d 3 q, we evaluate the doubles-to-singles correlation function 
as 




da^ 



(1) 
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C «/»(a a )- WW(2) ^^ ( 2) 
° l<7av ' 9relJ " WW ( qi ) W« (q 2 ) ' l j 

To allow an HBT analysis of C d//s , we have to make statistical assumptions (explained 
in the next section) that reduce the information contained in the two-particle quantities 
to that contained in the density matrix for a single particle source. The latter is conve- 
niently expressed in the Wigner representation as the (one-particle) phase space distribution 
gW(t, x,u, q), where uj is the energy of the meson. The connection to the correlation func- 
tion is given by Pratt's formula ||, which we will derive in the next section. The formula 
reads 

^Pratt/ \ = 1A- I d4xid4x2 9 (1) ( x l, 9av)ff (1) (x 2 , 9av) COS g re i(xi - X 2 ) 

l9av ' 9rolJ /d*x 1 <7W(x 1 ,ft)/d*a: 3 y(i)(x 2) fe) ' [) 

We use a notation here with four- vectors x = (t,x) and q = (u,q). To separate on-shell 
from off-shell four-momenta we use the notation q = (u(q),q) in the on-shell case, where 



uj(q) = J q 2 + m 2 is the energy of a meson with momentum q and mass m^. We have also 
defined the relative and the average four-momenta of the two mesons, q av = (q\ + 92)/2 and 

9rcl =91-92- 

To compare with theoretical expectations, it is common to model the evolution of the in- 
teracting hadronic system using classical transport equations P-Q] . With these simulations 
one can predict a classical source distribution function, g^ cl (t, x, \Jrri 2 + q 2 , q). Note that 
the energy is not an independent variable in the classical source. Thus in principle there is 
not enough information to use gW cl in Eq. ([[J), since the required g av in that equation need 
not have the energy on mass shell. However, in favorable cases the extrapolation off mass 
shell will not cause difficulties. 

We shall also find it useful to analyze the correlation data by dividing the two-particle 
distribution function by the unsymmetrized two-particle distribution. This is defined 

c i/u (q x WW( qi ,q 2 ) 

\q 1} q 2 ) 

where is the two-particle decay rate in the absence of Bose correlations. These corre- 
lations are not present for pairs of unlike mesons, so one could imagine measuring C l ^ n by 



e.g. comparing the correlation between pairs of ir + mesons to the correlation between 7r + 
and 7T~. We therefore shall call this the like-to-unlike correlation. However, although the 
Bose interference is absent from the unlike decay distribution, the interaction between unlike 
mesons induces stronger correlations than for charged mesons of like sign. Another approach 
to measure is to try to extract it from the correlated data by expressing the rate as a 
sum of two terms, one the uncorrelated product of two single-particle decays, and the other 
the correlation. This has been done in some experimental data analyses [|12|, [5|. We will see 
in this work that it is easier to satisfy the statistical assumptions for the like-to-unlike ratio 
than for the ratio of doubles to singles. 



II. DERIVATION OF THE PRATT FORMULA 

The basic formula for the correlation function in terms of the Wigner distribution of 
the source, Eq. (|3|), was derived by Pratt assuming an ensemble of classical currents. We 
will only use quantum mechanical amplitudes, as for example calculated from Feynman 
diagrams, as it is useful to see what assumptions are necessary to derive the formula from 
such quantum mechanical expressions. 

Our starting point is the symmetrized amplitude 

T s a ( qi , q 2 ) = T« ( qi , q 2 ) + (q 2 , q x ) (5) 

that describes the production of two identical mesons with momentum qi and q 2 in some 
reaction. Thus the events that contribute to the cross section o"o are those that produce 
just two mesons. The label a represents all final state variables that are not observed, such 
as the number of particles and their momenta. We separate these variables into three sets 
a = (ai|a2|a3) as follows. In the unsymmetrized amplitude ^( a i| a2 i a3 ) {Q11Q2) the label ot\ 
denotes the variables that directly influence the production of the meson with momentum 
qi, the label a 2 the ones that influence the meson with momentum q 2 , whereas 0:3 denotes 
all other variables in the final state. 
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Up to a possible overall sign, the unsymmetrized amplitude satisfies 

T (al\ a2 \a 3 ) (ft' ft) = T(a 2 \ ai \a 3 ) (ft> Ql) ■ (6) 

The differential decay rate of the initial state into two mesons is given by 

WP (q lt q 2 ) = \Y, \T S a (ft, ft)f (2tt) 4 5 4 (P a + q 1 + q 2 -P). (7) 

Here P denotes the energy-momentum vector of the initial state. P a is the energy-momentum 
carried away by all unobserved degrees of freedom. The sum over a represents summation 
over discrete degrees of freedom and integration over three-momenta, / d 3 p/(2n) 3 . The 
single-particle decay rate and W are given by integrals of over one and over both 
momenta, respectively. 

To obtain expressions for the correlation functions given in Eq. @ and Eq. (§) we first 
define the two-particle Wigner function 

9 g, (li , 9i;i2 , 92) ^ ?/ ^ e -, to ,, +Ml 

xTf ( gi + Pl /2, q 2 + pa/2) Tf * ( 9l - Pl /2, g 2 - p 2 /2) 

x(2n) A 5 4 (P a + qi + q 2 -P). (8) 

The notation Wp u for the normalization will be explained below. The inverse Wigner trans- 
form reads 

777- £ T ^ 2) (ft. ft) (ft' ft) (^ f ^ + E ft - ^ 

^-Pu a V i=l / 

= I A l( J 4 a; 2 e i[(91 -* +(92 - ? ^ 2l i 2) s a , *±*) . (9) 

Inserting (^) into (|7|) yields 
wf } (<7i,<7 2 ) 



J d 4 x 1 d 4 x 2 gp ] (xx, g x ; x 2 , g 2 ) + ^ (xi, g av ; x 2 , g av ) cos (g re i (x% - x 2 )) . 

(10) 

We can now express the like-to-unlike correlation function C^ u , Eq. (|]), in terms of two- 
particle Wigner functions if we identify the unlike particle production rate with the unsym- 

(2) 

metrized two-particle decay rate W P £ 



wPMx, 9 2 ) = E |^ 2) (ft, ») 2 (a*) 4 * 4 + ft + ft - P) 



d 4 xid 4 x 2 fi-p (xi, gi; ac 2 , 92) ■ 



If we now demand that g P be normalized for on-shell emission 



f d°q 1 f d°q 2 ( 2 ), _ ~M_ n 



(ii) 



we see that the normalization factor Wp M introduced in Eq. (^) is just the unsymmetrized 



total decay rate. We thus find for this correlation function 



C VU (<Z av ,<Z : 



rel 



)=1 + 



J d 4 x l d 4 x 2 g P ) (xi, g av ; x 2l q av ) cos(g re i(xi - x 2 )) 
/ d 4 xid A x 2 gp (xi, qi;x 2 ,q 2 ) 



(12) 



To get an expression for the correlation function in terms of one-particle quantities (c.f. 
Eq. (H)) the two-particle Wigner function must be factorized into one-particle functions. 
This demands a number of assumptions. In particular, we require that: 

(Al) the unsymmetrized T- matrix factorizes 



This condition can be fulfilled if the production process is dominated by tree graphs. 
However, this is not enough to factorize Eq. (H), because the 5- function still couples 
the aii and a 2 subsystems. So we also assume that: 

(A2) the summation over the degrees of freedom serves as a heat bath that effectively 
decouples the other two subsystems. 

Under these two conditions the two-particle Wigner function no longer depends explicitly on 
the energy and momentum of the initial state but depends on the statistical distribution of 
energy provided by the integration over 0:3. In effect, a microcanonical ensemble is replaced 
by a canonical ensemble characterized by a temperature in the center-of-mass frame 




(ft,ft)=T« ( ft )T« (q 2 )T c 



ibath 

a 3 



(13) 



]T(27r) 4 5 4 (P a + ft + q 2 - P) ~ exp((P - P° - P 2 ° - u x - u 2 )/T). 
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We defined Pf to be the energy for the degrees of freedom «j and uJi = u (q { ) = yqr? + ml 
is the energy of a meson with momentum q { . We label the corresponding Wigner functions 
9t, i.e., 

(1)/ v _ 1 f (1% 



#(^1) = ^E / + p 1 /2)T«*(g 1 - ft/2). (14) 

Then under (Al) and (A2) the two-particle Wigner function factorizes as 

g ?( Xl , qi ;x 2 , qi )= g ?( Xl , qi ) g ?(^ q2 ). (15) 

The requirements (Al) and (A2) correspond to the usually cited assumption that the 
meson field is chaotic, i.e., the mesons are created independently. Under these conditions 
we can now identify the like-to-unlike correlation function, Eq. (|]), with Pratt's expression, 
Eq. (HI). However, the two conditions above do not suffice to rewrite the doubles-to-singles 
correlation function, Eq. (0), in the same way. In addition we have to require that: 

(A3) the one-particle decay rate can be calculated neglecting the symmetrization of the 
amplitude. 

As this assumption is independent of (Al) and (A2) we can state it in terms of the micro- 
canonically calculated function as defined below Eq. (^) 



WjP fa) « Wj& = / ^% £ \U 2) (qi, g 2 )f (2vr) 4 tfW (P a + q 1 + q 2 - P) (16) 



A further consequence of (A3) is that the total decay rate Wp can be determined neglecting 
symmetrization as well, i.e, Wp ~ Wp u . It is usually assumed that this condition is fulfilled 
if typical momenta of the participating particles are large compared to the inverse size of 
the system, i.e., that the emission processes are sufficiently localized in the emission volume. 
This condition, however, can be misleading. We will come back to this point in Sec. [V]. 

Given assumption (A3) in addition to (Al) and (A2) the doubles-to-singles correlation 
function, Eq. (0), can be written using one-particle Wigner functions as in Pratt's expression, 
Eq. (§. 



III. A SOLVABLE MODEL 



In this section we will describe a simple model that we use to test how well the assump- 
tions (A1)-(A3) can be fulfilled in a quantum mechanical calculation of meson production 
by resonance decay. The model is based on a picture of a small source that emits heavy par- 
ticles. These undergo a cascade of decays to reach the final state. This picture is a plausible 
one for the decay of the concentrated high energy zone produced by nucleon-antinucleon 
annihilation or by e + -e~ annihilation at high energy. The mesonic source would have a 
spatial extension depending on the initial size of the state and the distance the resonances 
propagate. In heavy ion collisions, an additional mechanism extending the source is the 
rescattering of mesons in the later stages of the evolution. In our model, we will permit 
only the resonance propagation to contribute to the classical source size. The model is thus 
quite unrealistic for the heavy ion studies, but the criteria we develop have a more general 
validity. 

We will further simplify the resonance cascade process to the Feynman graphs shown in 
Fig. |l|. Here an excitation having definite energy and momentum decays into two resonances 
and N additional particles. Each resonance propagates and emits a meson into the final state. 
Only the two mesons are observed; the other particles will be integrated out in calculating 
the differential decay rates. Note that this graph satisfies the condition (Al) with the three 
sets of quantum numbers given by the momenta ki, k^, and the set pi. 

For computational convenience we assume that all unobserved particles have the same 
mass rrif, and ignore effects of symmetry among the unobserved particles. We further assume 
that we can treat the heavy particles nonrelativistically. Their on-shell energy is then given 
by e f (k) = k 2 /2m f . 

The intermediate resonances are assumed to be represented by simple resonance energy 
denominators of the form 




(17) 



S 



where e r (P) = P 2 /2m r , and m^m = m r — m/ is the difference between resonance mass 



and the mass of the unobserved particle in the final state. The quantity T r is the width of 
the resonance. In principle, the coupling constants should be such that the decay rates are 
consistent with the assumed widths, but since the absolute decay rates play no role we shall 
ignore this requirement and choose the coupling constant being one. 

The couplings at the vertices will be taken as point-like, since we are ignoring finite sizes 
of the initial source and of the particles. The meson-production vertices contribute factors of 



1/ yj2u(q) in the perturbation theory. For the numerical calculations, we choose parameters 
to correspond to the ir meson, emitted from the A resonance of the nucleon. Accordingly, 
= 139 MeV,m f = 939 MeV,m r = 1232 MeV, and T r = T A = 115 MeV. However, 
note that our assumption of constant form factors corresponds to an s-wave rather than the 
physical p-wave character of the A resonance. 

Finally, we will consider only one spatial dimension in the remainder of this work. We 
will switch to a notation with the space or momentum variables denoted with a Roman letter 
and the time or energy variable explicitly written out rather than subsumed in a relativistic 
notation. 

The differential decay rate associated with Fig. [I] is given by the following expression 





II -£ (2tt) 2 8{E -E 1 -E 2 -Y J e f (pi))S(P - P - P 2 - J» 



2 2tu(q 1 )2cu(q 2 ) J 2% 2% 



j=l i=l i=l 



x G r (E 1 ,P 1 )G r (E 2 ,P 2 ) 



2 



+ G r (E 1 + u(q 2 ) - u(qi), P\ + q 2 - qi)G r (E 2 + u(qt) - u(q 2 ), P 2 + qi - q 2 ) 




For notational convenience we have defined quantities 



Ei = e f (ki) +U}(qi) 



(19) 



and 



Pi = qi + h 



(20) 
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corresponding to the resonance energies and momenta, respectively. In the next section we 
will examine how this expression approaches Wt as the number of particles N gets large. 

IV. THERMAL LIMIT 

Our first application of the model will be to see how the thermal limit is approached 
as the number of spectator particles increases. We first define the spectator phase space 
volume by integrating the delta-functions in Eq. fll8|) over the spectator momenta pj, 

rt N) (E, P)=j dj ^0P L (2-) 2 S^E-J^f (ft)) 5 (p - £ ft) . (21) 
This may be evaluated in closed form to obtainf] 

JV-3 

^HE,P) = K{E-^-j 2 . (22) 

Here K is a constant that we calculated but, as the expression is complicated and will drop 
out later on, we do not show it here. Inserting this in Eq. ([iq), the two-particle probability 
distribution can be evaluated as a double integral. Before presenting the numerical results 
it is useful to see how the thermal limit arises. When N is large we may write the needed 
phase space function as 

qW( E - El -E 2 ,P-P 1 - P 2 ) « QW(E, P) exp(-N(E 1 + E 2 - (P 1 + p 2 )^—)/2E*), 

rrifJM 

where E* = E — P 2 /2rrifN is the center-of-mass energy. From here on we shall consider 
momenta in the frame where P = 0. The effective temperature associated with the phase 
space is then T = 2E/N, as expected from the equipartition theorem for N degrees of 
freedom. 

We now compare the doubles-to-singles correlation function evaluated explicitly from the 
finite- N two-particle probability function Eq. ( JlSf) , and the corresponding thermal limit, 



The iV-fold integration can only be done analytically for nonrelativistic particles. 
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w4 2) (?i,? 2 ) 

nW(£,P = 0) 1 rdhdk 2 



/-- exp( _ (£l+£2)/T) 



2 2cj(g 1 )2cj(g 2 ) •/ 2tt 2tt 

G r (E 1 ,P 1 )G r (E 2 ,P 2 ) 

+ G r (E t + u(q 2 ) - w(gi), Pi + q 2 - qi)G r (E 2 + cj(<?i) - u(q 2 ), P 2 + qi~ q 2 ) 



2 



(23) 



Fig. shows typical results for C d / S calculated with Eqs. ([18]) and (|23|) . The total energy 
was chosen to be E = 1.5m n N corresponding to a temperature of T = 3m„ which is 
the temperature used in Eq. fl2"3"|). The correlation function is shown as a function of the 
momentum difference q re \ = q\ — q 2 for fixed average momentum q av = 0. It may be seen that 
the finite- N correlation function approaches the result of the thermal limit. So we conclude 
that as few as 50 participating particles justify the assumption that the correlation function 
can be calculated using a canonical ensemble instead of a microcanonical one. We expect 
that in three dimensions the same heat bath would require of the order of 20 particles. 



V. THE SINGLES DISTORTION 

It is immediately apparent from Fig. |2| that the correlation function looks quite different 
from the expected form, which should start at 2 at zero q re \ and fall monotonically to 1 
at large q re \. We can trace this behavior back to a violation of condition (A3). In Fig. |3| 
we show Wt (full line) calculated by integrating over the symmetrized as compared 
to an integral over the unsymmetrized (dashed line) (cf. Eqs. (|7|) and (|TT)D). The 

two curves clearly deviate strongly for small values of q. It is obvious that also the total 
particle production probability can not be calculated neglecting symmetrization. These facts 
account for the wrong asymptotic values of the correlation function. To make things worse 
the correlation function has a pronounced dip for intermediate values of q re \. That means 
that even the introduction of an overall normalization factor and a coherence parameter A 
does not help to extract the desired information about the source from C d//s . 

The validity of (A3) can be determined empirically . The single particle distribution 
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is evaluated by integrating over one of the two momenta in the two-particle distribution 
function Eq. ( |I~0D or (p3|). The contribution of the interference term is of the order of the 
width of the correlation function Aq c divided by a typical momentum difference 
< (Qi ~ Q2) 2 > 1 ^ 2 - This should be small for (A3) to hold. From Fig. ^| we have Aq c « 
and Fig. [3] gives < (qi — q 2 ) 2 >^ 2 ~ Sm^, yielding 

, Aq \ 2 m « 25%. (24) 
< (9i - > 1/2 

This is intolerably large to ignore in the HBT analysis. One might anticipate that in three 
dimensions the relevant expression would be the cube of Eq. fl24|) , which would make (A3) 
a rather good assumption. However, analysis of experimental data shows distortions of the 



order of 10 % H . Also in a Fritiof Monte Carlo simulation of NA35 data [|1^] an influence 
of the symmetrization of the singles spectrum of 10-20 % was found. Coming back to our 
model, we find that by decreasing the width of the resonance down to T r = Ta/10 we actually 
decrease Aq c so much that (A3) becomes valid. In Fig. |] we show the single particle decay 
rate again for the symmetrized (full line) and unsymmetrized (dashed line) two particle 
amplitude. Compared to Fig. ^ the deviation of the two curves is very small. Indeed the 
calculation of the doubles to singles correlation function yields a satisfactory result, shown 
in Fig. [| (full line). The small deviation from two at low g re i and one for high q re \ is due 
to the residual small difference between the single particle decay rate in the symmetrized 
and the unsymmetrized case respectively as well as numerical inaccuracy in evaluating the 
two-particle decay rate for such a small resonance width. In Fig. |5| we show also the like-to- 
unlike correlation function C 1//u for the same parameters (dashed line). As expected for the 
case of a small singles distortion we find a good agreement of the two correlation functions. 

VI. SOURCE SIZES 

In the remainder of this work we shall concentrate on the measure C l ^ n to avoid the 
problem of the distortion in the singles spectrum. In Fig. we show the thermal model 
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q\/m £ or p arame ters corresponding to the physical A resonance and for several values of 
the temperature. The correlation function is seen to be well-behaved, and one can try to 
extract source size parameters from the width of the peak. Notice that the peak becomes 
the narrower as the temperature increases. This seems to indicate that the emitting system 
is larger for higher temperatures. We will discuss this point in detail in Sec. |VTB| . We now 
address the question of whether a classical source size can be extracted from the peak width. 



A. Classical limit 

We first derive the classical limit of the one-particle Wigner function for our model. The 
Wigner transform of the quantum amplitude has the form 

= f *-K+e /(fc ))/T f ^^ e -iM-«) (2w (gi) 2u (ft)r va 
J 2n J 2tt 2ir 

xG r (6f (k) + W av + k + ?l) G* (6f (k) + UJ m - Wrel/2, k + q 2 ) ■ (25) 

The integration over uj tc \ can be performed by an ordinary contour integration. The 
contribution from the two poles is 

= &(t)J * e -K +e ^)/ T J ^e^* (2u ( gi ) 2u (q 2 )y 1/2 
2ie~ rt 

X E +E (exp (-2iE + t) - exp {2iE_t)) , (26) 

where 

E± = e r (k + q av ± <7rei/2) + m diS - e f (k) - u av . (27) 

Up to this point the calculation is still exact, i.e., Eq. ( p6|) still contains the full quantum 
mechanical information. We now derive a classical source function by approximations on 
the above expression. We first expand E± to first order in the momentum difference q re \ 

de r (p) 



E±& AE± 



dp 
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qrei/2, (28) 



where 

AE = e r (k + g av ) + m diS - e f (k) - u av . (29) 

The derivative de r (p) /dp is the classical velocity v r (p) of the resonance. Thus we may write 

n {1)cl (t r n a ) - 9 ® r~ r ^ f — p -K+e/(fc)) /T f ^el i( grj3! - q to r (fc+g„)) sin (gAgg 
w o OT J 27T J 2tt AE 



(30) 



The integral over g av gives a 5-function and the last factor approaches the 5-function 
27r5(A£'). Therefore the classical approximation amounts to demanding that the resonance 
propagates on-shell. Evaluating the 5-functions and dropping the subscript "av" we obtain: 

© (t) -F r t X „. flJ , „\ A „-(k* f +q) 2 /2m r T 1 



# cl (t, x, Wj g) = ^±e~™ u r (k) + q) t) e~^) W 

^ (?) V V 7 1 V r (k) + q)-V f (k)) 

(31) 

In this equation, kj and kj are the two roots of the condition AE = 0, which resulted from 
the evaluation of the k- integral over the delta-function 6(AE). These roots are solutions of 
a quadratic equation and are given by 



k y = HlL q ± !^fp: q 2 + 2mf m r ( 1 - — ) . (32) 



m Am - y m5 ifr - \ m diS/ 

The k-integral gives the density-of-states factor at the end of equation (^l|). Note that the 
various terms in this Wigner function (except the density-of-states factors) have obvious 
classical interpretations. The initial G function specifies that the source starts at t = 0, and 
its exponential decay is given by the second factor. The resonance propagates classically with 
connection between its position, velocity, and the time given by the argument of the delta- 
function. The probability to make the resonance is given by the Boltzmann factor following. 
With this one-particle Wigner function in the classical approximation we can evaluate the 
classical correlation function C cl by inserting , Eq. (|3T|), into Pratt's expression, Eq. 



Note that our classical g? depends on the variables q and uj independently. Fig. |7| shows 
the correlation as a function of q re \ with average momentum fixed at g av = and temperature 
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T = m,,-, indicated by the solid curve. The correlation function behaves as expected for small 
momentum differences. For large momenta, however, C cl develops a singularity. This can 
be traced back to the occurrence of u av in the argument of the Wigner functions in the 
numerator of Eq. (|3]). Under sufficiently extreme conditions, this energy goes off-shell so far 
that the classical condition AE = no longer has a solution. If u av is replaced by an on-shell 
energy, namely by u(q av ), the problem does not arise. The correlation function calculated 
with this on-shell prescription is shown as the dashed curve. This function is very similar to 
the original one for small q ie \. In contrast to the singular behavior of the original correlation 
function the dotted curve behaves smoothly also for high q re \. Since this problem arises at 
high momentum, and the source size question depends on the low momentum behavior only, 
we can ignore the differences for the remainder of this work. 



B. Contributions to the source size 



A convenient way to characterize the information contained in the correlation function 
is to take the second derivative with respect to the momentum difference at q re \ = 0. This 
gives an effective mean square source size defined as 

1 d 2 C 



< X 2 > Cff E 



(33) 



<3rol=0 



2 dq 2 el 

We determine this from the single-particle Wigner source function taking the derivative of 
Eq. (|). This yields 



< x 2 > cS 



^ J dtdx (x - v n (q w ) tf g- J dtdx (x - v w (g av ) t) g 
1 d 2 G 1 d 2 u{q) 



' 1 dG 
2Gdq av 



AG dql AG dq 2 



9=<?a 



d_ 



dtdx g{t, s^g^U^) 



(34) 



In writing this we used the abbreviations g = g^ (t, x, u>(q av ), g av ) and G = 
J dtdxg^ (t, x, ui(q av ), q&v)- The quantity v n is the velocity of the pion v n (q) = du(q)/dq. 
We will refer to the five contributions of the r.h.s. of Eq. fl3"l| ) by < x 2 > 7 , . . . , < x 2 > v re- 
spectively. < x 2 > I and < x 2 > H correspond to the space-time distribution of the emission 
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points of the mesons. This is the information we would like to extract. The other terms 
give corrections which are due to the energy and momentum dependence of the source. 

To simplify the discussion we concentrate on the behavior at g av = 0. Then < x 2 > !I 
and < x 2 > IH in Eq. fl53|) vanish by symmetry. The first term is easily evaluated with the 
classical source g4 (Eq. ([H])). The result is 

\2 (l)cl 



1 f 

<x 2 >[ x = — J dtdx (x - v n (q av ) t) 2 gf' 

(35) 



v r / 

where k® is the classical momentum of the resonance required to emit a meson at q = 0. It 
is given by (compare Eq. ([32])) 



The result Eq. ( pop is of course exactly what one would expect for a source moving with 
velocity k®/m r and decaying at a rate T r . 

Note also that the source size is independent of the temperature. This is counterintuitive 
at the first glance. However, in the classical approximation the momentum of the resonance 
is kinematically fixed by the momentum of the emitted pion. Therefore the distance the 
resonance travels before it emits the pion is determined only by the kinematics and lifetime. 

C. Comparison of classical and quantum results 

We now extract the source size from the quantum mechanical thermal source. The 
correlation function at q av = has the form 

(g av = 0, q iel ) = 1 + IM^-fteOI 2 (37) 

where 

X T ( qi ,q 2 ) = J dk e -M*)-MW2))/T 

xG r (e/(A;) + co( qi /2), k + q x /2) G* r {e f {k) + uj{q 2 /2), k + q 2 /2). (38) 
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It is straightforward to evaluate X by contour integration. The result in the T — > oo limit is 
particularly simple, yielding for the correlation function 

Ci u (g av = 0, g rel ) = 1 + -2 , (39) 

(l + (g re i/go) J 

where we have defined 



g = Im. ^q 2 el + (m diff - U (g rel /2) - ,T P /2). (40) 

V m / m f 

The source size may be evaluated from Eq. (^). In the limit T r -C rridis — win, it reduces 
exactly to the classical formula, Eq. fl35|), 



< x 2 >f «< x 2 >i, . (41) 

At finite temperature, the expression for C 1//u is rather unwieldy and we quote only the 
source size, obtained either directly from Eq. ( |3"B| ) or with Eq. and the quantum Wigner 
function Eq. (^5J). The quantum mean square radius < x 2 is found to be 



< x 2 >! ff - 



7 (v r -v n y 



V r - V m 



\AE + iT r /2 i 



(42) 



q \\(AE) 2 + (r r /2) 2 / / 

Here we used AE = e r (A;+g a v)+m ( jifr— £/(&) — uj(q av ) and the thermal average ((/)) weighted 
by the resonance distribution is defined as 

/ dk e~ e f^/ T f ( (AE) 2 + (T r /2) 2 ) ^ 

({f))= ~ ; 2 2 L ■ (43) 

fdk e-/( fc )/ T [(AE) 2 + (r r /2) 2 ) 

In Fig. ^| we show the mean square source size for a width of the resonance T r = Ta 
and g av = 0. Note that for vanishing average momentum only the first term on the r.h.s of 



Eq. (H2T) contributes. The quantum mean square radius < x 2 >^ ff is shown as a function of 



9 

temperature, indicated by the solid line. This is compared to the temperature-independent 
classical size (x 2 )^ indicated by the dashed line. It might be hoped that the additional 
terms in Eq. (|34D would improve the classical description, but this is not the case. The full 
classical source size < x 2 >lf , shown as the dot-dashed line, in fact deviates in the opposite 
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direction from the quantum result. The disagreement is serious, since we are interested in 
temperatures of the order of T = m n and lower, for physical sources. 

Clearly, the off-shell propagation is important at physically realizable temperatures. Even 



small width, T r = Ta/10. We see that even in this case the classical source is unreliable 
for temperatures below 2m n . The problem is severe because in the classical calculation a 
zero-momentum pion is not produced from a low momentum resonance as is favored in the 
quantum mechanical calculation. However, if we choose conditions such that the resonance 
region is not suppressed by the Boltzmann factor, the classical source is accurate down to 
much lower temperatures. 

Even though our model calculation shows that quantum effects are important for the 
interferometry among mesons produced by resonances, our result Eq. ( [42|) suggests a way to 
improve classical modeling. When particles propagate far off-shell, their lifetime is controlled 
by (AE)" 1 , according to the uncertainty principle, rather than by T" 1 . Eq. ( f£2| ) could be 
recovered from a classical simulation if in Eq. ( |55| ) the resonance lifetime T" 1 is replaced by 
an effective lifetime ~ ((r r /2) 2 + (AE) 2 )^ 1 ^ 2 , with a suitable averaging over an ensemble of 
resonances with different A£"s. 



We have studied a resonance decay model to test the statistical assumptions behind the 
HBT analysis of meson correlations. One set of assumptions gives rise to a reduction of 
the unsymmetrized two-particle source function to a product of two single-particle thermal 
sources. We find that it is not too difficult to realize conditions that substantially satisfy this 
statistical assumption. The usual HBT analysis requires, however, that the single-particle 
distribution be calculable from an unsymmetrized two-particle source, and this condition is 
much more difficult to fulfill. 




VII. SUMMARY 
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For resonance decays, we can ask whether the correlation reflects the spatial propagation 
of the resonances before they decay. We find that the apparent source size extracted from the 
quantum mechanically calculated correlation function deviates strongly from the classical 
mean square size even for resonance parameters and temperatures that suggest a classical 
behavior. The off-shell propagation of resonances is important and needs to be taken account 
of in the classical propagation if extracted source sizes are to make sense. 
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FIGURES 

FIG. 1. Graphical representation of the symmetrized amplitude for the production of two 
mesons in our model. 

FIG. 2. Correlation function C d / S for different numbers of particles and in the thermal limit. 
The energy is E = 1.5Nm n which corresponds in the thermal limit to a temperature of T = < &m- K . 
The average momentum of the meson pair is zero for all cases. 

FIG. 3. Single particle decay rate W^(q) calculated with the symmetrized (full line) and 
unsymmetrized (dashed line) two-particle amplitude, for T = 3m n and T r = 

FIG. 4. Same as Fig. ^, but for a small resonance width T r = Fa/10. 

FIG. 5. The correlation functions C d / S (full line) and C l / n (dashed line) for a resonance width 
T r = Ta/10 and T = Sm^- 

FIG. 6. The correlation function C l ^ n for different temperatures T = 0.5m n (full line), LOm^ 
(dashed line), and fy.Om^ (dot-dashed line). The width of the resonance is T r = Ta- 

FIG. 7. Classical correlation function with different treatments of the energy variable, for 
T = 1 m-x and r r = Ta- 

FIG. 8. The mean square source size as a function of the temperature in various treatments. 
The horizontal dashed line is the classical size < x 2 >^, Eq. (^). The dot-dashed curve includes 
corrections to the size of the classical source from Eq. (^). The solid line gives the result from the 
quantum mechanical calculation, Eq. (|42]). The width is T r = Ta- 

FIG. 9. Same as Fig. || but for the reduced width T r = Ta/10. 
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Figure 1; G.F. Bertsch 
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